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Abstract 

Cell  temperature  control  plays  a  crucial  role  in  SOFC  operation.  In  order  to  design  effective  temperature  control  strategies  by  model-based  control 
methods,  a  dynamic  temperature  model  of  an  SOFC  is  presented  in  this  paper  using  least  squares  support  vector  machines  (LS-SVMs).  The  nonlinear 
temperature  dynamics  of  the  SOFC  is  represented  by  a  nonlinear  autoregressive  with  exogenous  inputs  (NARXs)  model  that  is  implemented  using 
an  LS-SVM  regression  model.  Issues  concerning  the  development  of  the  LS-SVM  temperature  model  are  discussed  in  detail,  including  variable 
selection,  training  set  construction  and  tuning  of  the  LS-SVM  parameters  (usually  referred  to  as  hyperparameters).  Comprehensive  validation  tests 
demonstrate  that  the  developed  LS-SVM  model  is  sufficiently  accurate  to  be  used  independently  from  the  SOFC  process,  emulating  its  temperature 
response  from  the  only  process  input  information  over  a  relatively  wide  operating  range.  The  powerful  ability  of  the  LS-SVM  temperature  model 
benefits  from  the  approaches  of  constructing  the  training  set  and  tuning  hyperparameters  automatically  by  the  genetic  algorithm  (GA),  besides  the 
modeling  method  itself.  The  proposed  LS-SVM  temperature  model  can  be  conveniently  employed  to  design  temperature  control  strategies  of  the 
SOFC. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  solid  oxide  fuel  cell  (SOFC)  is  expected  to  be  one  of 
the  most  promising  devices  for  energy  conversion  because  of 
its  high  efficiency,  low  pollutant  emissions  and  flexible  fueling 
strategies.  Possible  applications  of  SOFCs  range  from  vehicular 
auxiliary  power  units  to  stationary  power  plants  [  1  ] .  Cell  temper¬ 
ature  control  plays  a  crucial  role  in  SOFC  operation  [2,3] .  This  is 
because  the  cell  temperature  has  a  significant  effect  on  the  elec¬ 
trical  characteristics  of  an  SOFC,  e.g.  the  cell  voltage  and  the 
current  density  distribution.  Moreover,  the  variation  of  cell  tem¬ 
perature  may  cause  degradation,  even  damage  to  the  fuel  cell  due 
to  thermal  fatigue  or  thermal  cracking  of  the  cell  material.  Thus, 
the  average  cell  temperature  should  be  carefully  controlled.  In 
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order  to  design  effective  temperature  control  strategies  by  using 
model-based  control  methods,  simple  and  accurate  dynamic 
temperature  models  of  SOFCs  are  highly  desired. 

In  recent  years,  many  dynamic  physical  models  of  SOFCs 
based  on  conservation  laws  and  electrochemical  principles  have 
been  reported  [3-6] .  Although  these  models  are  useful  for  the 
analysis  and  optimization  of  SOFCs,  in  general  they  are  too 
complex  for  applying  model-based  control  methods.  The  com¬ 
plicated  phenomena  associated  with  SOFCs  easily  drive  such 
models  to  high  complexity  [7] .  Developing  simpler  data-driven 
models  based  on  the  system  identification  methodology  is  an 
important  way  to  control  complex  plants  [8,9].  Following  this 
line,  many  data-driven  modeling  approaches,  such  as  classi¬ 
cal  system  identification,  artificial  neural  networks  (ANNs)  and 
fuzzy  logic,  can  be  applied  to  dynamic  temperature  modeling 
of  SOFCs.  The  requirements  for  developing  such  data-driven 
temperature  models  of  SOFCs  for  control  purpose  can  be  sum¬ 
marized  as: 
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Nomenclature 

A 

cell  active  area  (m2) 

Ar 

reforming  reaction  surface  area  (m2) 

^sh 

anode  section  area  (m2) 

AR 

air  ratio 

b 

bias  term 

bo-b?, 

fitting  function  coefficients  of  specific  heat  capac¬ 
ity 

cp 

specific  heat  capacity  of  gas  species 
(Jmol-1K-1) 

c*p 

specific  heat  capacity  of  solid  part  (Jkg-1  K-1) 

ek 

prediction  error  of  the  Mi  data  point 

E* 

activation  energy  (J  mol-1) 

F 

Faraday’s  constant  (=96,485  C  mol-1) 

hi 

molar  enthalpy  (J  mol-1) 

AH 

total  change  of  molar  enthalpy  for  chemical  reac¬ 
tion  (Jmol-1) 

J 

average  current  density  (A  m-2) 

kp 

coefficient  for  unit  conversion  (=105  Pa  bar-1) 

kr 

reforming  reaction  constant 

(=4274  mol  s-1  m-2  bar-1) 

ksh 

WGS  reaction  constant  (=1.2  x  104molm-3  s-1) 

K 

kernel  function 

equilibrium  constant  of  WGS  reaction 

M8 

mass  of  solid  part  (kg) 

number  of  transferred  electrons 

m 

molar  number  of  species  i  within  SOFC  (mol) 

ny ,  rid ,  n 

u  maximum  lags 

N 

molar  flow  rate  (mol  s-1),  number  of  data  points 

P 

pressure  (bar) 

P  dc 

direct  current  electric  power  (W) 

ri 

rate  of  reaction  j  (mol  s-1) 

R 

universal  gas  constant  (=8.314  Jmol-1  K-1) 

^ohm 

total  cell  resistance  (£2  m2) 

T 

temperature  (K) 

To 

initial  temperature  (K) 

U 

voltage  (V) 

Uf 

fuel  utilization 

V 

volume  (m3) 

w 

feature  vector 

X 

input  vector 

Xi 

molar  fraction  of  species  i 

y 

output  value 

Greek  symbols 

Oik 

Fagrange  multiplier 

fian>  Pea 

pre-exponential  factors  of  anode  and  cathode 
(Am-2) 

11 

overpotential  (V) 

<P 

nonlinear  mapping 

y 

regularization  parameter 

(7 

kernel  parameter 

^an?  ^ca 

,  (7 eiec  electrical  conductivity  of  anode,  cathode 
and  electrolyte  m-1) 

r 

thickness  (m) 

§  number  of  chemical  reactions 

f  number  of  gases 

Subscripts 

act  activation 

an  anode 

ca  cathode 

cone  concentration 

elec  electrolyte 

i  gas  species  i 

j  the  jth  reaction 

k  sampling  instant  k ,  the  Mi  data  point 

ohm  ohmic 

r  steam  reforming  reaction 

sh  water-gas  shift  reaction 

Superscripts 
in  fuel  cell  inlet 

out  fuel  cell  outlet 

ref  reference  condition 

s  solid  part  (including  the  PEN  and  interconnect) 


(1)  The  model  should  be  able  to  represent  the  temperature 
dynamics  of  an  SOFC  accurately  over  a  relatively  wide 
operating  range,  certainly  including  the  design  point. 

(2)  The  model  can  be  conveniently  integrated  into  the  whole 
control  scheme  of  an  SOFC. 

(3)  The  model’s  performance  should  be  fully  tested  over  the 
expected  operating  range. 

(4)  The  modeling  process  should  be  as  simple  as  possible. 

Compared  with  the  substantial  progress  on  physical  model¬ 
ing,  however,  data-driven  dynamic  temperature  models  of  fuel 
cells  (particularly  of  SOFCs),  which  meet  all  above  require¬ 
ments,  are  still  rare. 

Shen  et  al.  [10]  presented  an  application  of  radial  basis  func¬ 
tions  (RBFs)  neural  networks  to  develop  a  dynamic  temperature 
model  of  a  molten  carbonate  fuel  cell  (MCFC)  stack.  Although 
satisfactory  degree  of  precision  is  obtained,  the  MCFC  RBF-NN 
model  has  several  weaknesses.  First,  the  topology  and  train¬ 
ing  strategy  of  RBF-NNs  are  determined  by  experience,  which 
weakens  the  model’s  objectivity.  Second,  RBF-NNs  get  easily 
stuck  in  local  extrema.  Third,  the  effect  of  load  disturbance  on 
cell  temperature  was  left  out  of  consideration,  which  impairs  the 
model’s  generality.  Jurado  [11,12]  used  different  Hammerstein 
models  to  represent  the  dynamic  characteristics  of  an  SOFC. 
As  the  Hammerstein  model,  which  uses  a  static  nonlinear  block 
followed  by  a  dynamic  linear  block  to  represent  nonlinear  sys¬ 
tems,  is  a  special  type  of  nonlinear  model,  these  SOFC  models 
still  need  more  tests  over  a  wide  operating  range.  Yang  et  al. 
[13]  reported  a  dynamic  T-S  fuzzy  model  of  a  MCFC,  which 
consists  of  a  voltage  and  a  temperature  submodel.  However,  the 
temperature  submodel  can  only  be  valid  within  a  relatively  nar¬ 
row  temperature  range.  In  addition,  the  identification  procedure 


Y.-W.  Kang  et  al.  /  Journal  of  Power  Sources  179  (2008)  683-692 


685 


of  the  T-S  model,  which  contains  an  antecedent  identification 
and  a  consequent  identification,  is  complicated. 

As  an  important  breakthrough  in  the  area  of  machine  learn¬ 
ing  in  the  last  decade,  the  support  vector  machine  (SVM)  has 
proven  to  be  a  powerful  tool  for  solving  classification  and  regres¬ 
sion  problems.  However,  the  training  of  SVM  needs  to  solve  a 
quadratic  programming  (QP)  problem,  which  is  time-consuming 
in  the  case  of  large  data  set.  The  least  squares  support  vector 
machine  (LS-SVM),  which  achieves  the  solution  by  solving  a 
set  of  linear  equations,  is  a  least  squares  version  of  standard 
SVM  [14].  Compared  with  SVM,  it  can  significantly  reduce  the 
training  time  at  the  price  of  just  a  small  precision  loss  [14] .  More¬ 
over,  it  maintains  SVM’s  advantages  of  superior  generalization 
ability  and  finding  global  solution.  In  recent  years,  the  LS-SVM 
method  has  aroused  interests  of  many  researchers  [15-17],  and 
has  already  been  used  to  model  the  static  voltage  characteristics 
ofanSOFC  [18]. 

In  this  paper,  a  dynamic  temperature  model  of  an  SOFC, 
which  aims  to  satisfy  all  the  above  modeling  requirements,  is 
built  using  the  LS-SVM  regression  method.  The  nonlinear  tem¬ 
perature  dynamics  of  the  SOFC  is  represented  by  a  nonlinear 
autoregressive  with  exogenous  inputs  (NARXs)  model  [9]  that 
is  implemented  using  an  LS-SVM  regression  model.  In  this 
study,  the  training  data  is  generated  with  a  physical  model  of 
an  SOFC.  A  training  set  covering  the  entire  expected  operating 
range  is  constructed.  To  improve  the  resulting  model’s  general¬ 
ization  ability,  an  innovative  optimization  algorithm  known  as 
the  genetic  algorithm  (GA)  [19,20]  are  adopted  to  automatically 
tune  two  parameters  in  LS-SVM  design  (usually  referred  to  as 
hyperparameters).  In  addition,  the  LS-SVM  temperature  model 
is  fully  tested  over  the  operating  range  by  comparison  with  the 
physical  model. 

The  rest  of  this  paper  is  organized  as  follows.  In  Section  2,  the 
SOFC  process  is  briefly  described.  In  Section  3,  the  development 
of  the  dynamic  temperature  model  using  LS-SVMs  is  presented 
in  detail.  Modeling  results  and  model  validation  are  presented 
in  Section  4.  Finally,  conclusions  are  drawn  in  Section  5. 


Fuel 


|  Interconnect 

CH4  +  H,0  — >  CO  +  3H2 

co+h2o<->co2+h2 


Depleted 

Fuel 


Fig.  1.  Schematic  diagram  of  SOFC  process. 


Table  1 

Definitions  of  fuel  utilization  and  air  ratio 


Fuel  utilization 


Air  ratio 


Uf  = 


JA/2F 


AR  = 


( —I—  Y^  —1—  Y  ^ 

^ach4  ^  ah2  ^  V 

ill  yyin 


W^uel 


02  vair 

JA/AF 


can  be  described  as: 

Anode:  H2  +  02~  H20  +  2e~,  (la) 

Cathode  :  \C>2  +  2e“  O2-.  (lb) 

For  a  direct  internal  reforming  SOFC  (DIR-SOFC)  that  operates 
on  methane,  the  internal  reforming  process  also  takes  place  in  the 
anode  channel  to  produce  hydrogen,  and  is  usually  represented 
by  the  following  steam  reforming  reaction  and  water-gas  shift 
(WGS)  reaction: 

Steam  reforming  :  CH4  +  H2O  ->  CO  +  3H2,  (2) 

WGS  :  CO  +  H20  C02  +  H2.  (3) 


2.  Description  of  SOFC  process 

The  SOFC  can  be  viewed  as  a  electrochemical  reactor  oper¬ 
ating  at  very  high  temperatures  (600-1000  °C),  which  produces 
electricity  and  heat  directly  from  the  electrochemical  combina¬ 
tion  of  a  gaseous  fuel  (hydrogen  or  hydrocarbons,  e.g.  methane) 
with  an  oxidant  (typically  air).  An  SOFC  consists  of  an  intercon¬ 
nect  structure  and  a  tri-layer  structure  composed  of  two  porous 
ceramic  electrodes,  anode  and  cathode,  separated  by  a  dense 
ceramic  electrolyte  (often  referred  to  as  the  PEN).  The  operat¬ 
ing  principle  of  an  SOFC  process  may  be  described  in  Fig.  1. 
Fuel  and  air  are  fed  into  the  anode  and  cathode  channels,  respec¬ 
tively.  At  the  cathode,  oxygen  is  consumed  to  form  oxygen  ions 
by  obtaining  electrons.  The  oxygen  ions  migrate  through  the 
ion-conducting  electrolyte  to  the  anode  where  they  combine 
with  hydrogen,  producing  water  vapor  and  releasing  electrons 
to  the  external  circuit.  The  released  electrons  are  transferred  to 
the  cathode  via  an  external  load,  where  they  perform  electrical 
work.  The  electrochemical  reactions  for  the  anode  and  cathode 


The  cell  voltage  is  related  to  the  concentration  of  electro¬ 
chemical  species  by  the  well-known  Nernst  equation.  In  order 
to  provide  a  voltage  greater  than  zero,  not  all  the  reactants  in 
the  cell  can  be  consumed  [21].  Thus,  fuel  and  air  in  excess 
of  the  required  amount  of  electrochemical  reactions  should  be 
delivered.  In  addition,  since  the  electrochemical  combination  of 
hydrogen  and  oxygen  is  a  highly  exothermic  reaction,  excess 
air  is  also  needed  for  cell  cooling.  Two  performance  factors,  i.e. 
fuel  utilization  and  air  ratio  (see  Table  1)  are  usually  used  to 
define  the  amount  of  utilized  fuel  and  excess  air,  respectively. 

3.  Development  of  the  dynamic  temperature  model  of 
SOFC  using  LS-SVMs 

In  this  section,  the  LS-SVM  regression  method  is  first 
described.  Then  issues  concerning  the  development  of  the 
dynamic  temperature  model  are  discussed  in  detail,  including 
variable  selection,  training  set  construction  and  hyperparameter 
tuning. 
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3.1.  LS-SVM for  regression 


3.2.  Variable  selection  of  the  dynamic  temperature  model 


The  nonlinear  modeling  problem  based  on  measurable  data 
essentially  is  a  nonlinear  regression  problem.  The  LS-SVM 
regression  model  [14,22]  employs  the  following  linear  represen¬ 
tation  in  a  so-called  higher  dimensional  feature  space  H  c  Rnh : 

yO)  =  wT<p(x)  +  b ,  (4) 

where  x  e  Rn  is  an  input  vector,  y  e  R  is  an  output  value,  w  e  H 
is  a  feature  vector.  ^(-)  :  Rn  Rnh  is  a  nonlinear  mapping  from 
the  input  space  to  the  feature  space,  by  which  a  nonlinear  regres¬ 
sion  in  the  input  space  is  converted  to  a  linear  regression  in  the 
feature  space. 

Given  a  training  set  {xk,  yk)k=v  the  LS-SVM  method  con¬ 
siders  the  regression  problem  as  the  following  optimization 
problem: 


1  T 

min  J(w,  e)  =  -w  w 

w,b,e  2 


subject  to  equality  constraints 

yk  =  wT(p(xk)  +  b  + ek,  k  =  l, . . . ,  N, 


(5) 


(6) 


where  y  is  the  regularization  parameter,  ek  is  the  prediction  error 
for  data  point  k. 

By  applying  the  Lagrangian  multiplier  method,  the  solution 
to  above  optimization  problem  is  obtained  as  the  following  linear 
equations: 


0 

1 


iT 

Q  +  y~lI 


(7) 


where  a  =  [«i ,  . . «:v]T  is  the  Lagrange  multiplier  vector, 
y  =  [yu  ;w]T,  1  =  [1,  1  ]TeRN,  and  flisan]Vx/V 

kernel  matrix.  By  using  the  kernel  trick  [22],  one  obtains 
Qkl  -  <P(xk)T =  K(Xk,  xi ),  k,  1=1,  . . .,  N.  And  the  resulting 
LS-SVM  regression  model  becomes 


N 

y(x)  =  yfakK(x,  xk)  +  b, 
k= 1 


(8) 


where  c^,  b  are  the  solution  to  Eq.  (7). 

Note  that  the  dot  product  <p(-)T(p(-)  in  the  feature  space  is 
replaced  by  a  prechosen  kernel  function  K(  • ,  • )  due  to  the  employ¬ 
ment  of  the  kernel  trick.  Thus,  there  is  no  need  to  construct  the 
feature  vector  w  or  to  know  the  nonlinear  mapping  </?(•)  explic¬ 
itly.  Given  a  training  set,  the  training  of  an  LS-SVM  is  equal  to 
solving  a  set  of  linear  equations  as  Eq.  (7).  This  greatly  simplifies 
the  regression  problem. 

The  chosen  kernel  function  must  satisfy  the  Mercer’s  condi¬ 
tion.  The  linear,  polynomial  and  radial  basis  functions  kernel  are 
the  commonly  used  kernel  functions.  In  this  paper,  the  popular 
RBF  kernel  is  used  and  this  is  of  the  form 


In  order  that  the  resulting  dynamic  temperature  model  can  be 
integrated  into  the  whole  control  scheme  of  an  SOFC,  its  vari¬ 
ables  should  be  selected  first  according  to  the  employed  control 
scheme. 

The  aim  of  controlling  an  SOFC  is  to  deliver  specified  electric 
power  to  the  load  and  avoid  any  operating  conditions  that  may 
cause  degradation  or  damage  to  the  fuel  cell.  Due  to  the  limi¬ 
tations  of  components’  mechanical  properties  and  the  intrinsic 
nature  of  establishing  the  cell  voltage  mentioned  above,  the  aver¬ 
age  cell  temperature  and  the  supply  of  fuel  and  air  are  the  two 
most  important  aspects  that  should  be  controlled.  As  the  average 
cell  temperature  cannot  be  measured  conveniently,  in  practice  a 
measurable  variable,  i.e.  the  outlet  fuel  temperature  is  controlled 
instead  [2,3].  The  supply  of  fuel  and  air  can  be  reflected  by  the 
fuel  utilization  and  air  ratio,  respectively.  Hence  the  common 
control  objectives  of  an  SOFC  include:  maintaining  the  fuel  uti¬ 
lization  as  constant  as  possible;  ensuring  the  air  ratio  in  a  certain 
range;  guaranteeing  the  outlet  fuel  temperature  as  constant  as 
possible  [2]. 

These  control  objectives  can  be  achieved  using  the  control 
scheme  illustrated  in  Fig.  2.  Since  the  average  current  density 
is  decided  by  uncontrollable  load,  it  is  taken  as  a  disturbance 
variable  in  the  control  scheme.  Fuel  and  air  flow  rates  are  two 
basic  operating  degrees  of  freedom.  By  measuring  the  average 
current  density  as  a  feedforward  signal,  the  fuel  flow  rate  cal¬ 
culated  according  to  the  definition  of  fuel  utilization  can  keep 
the  fuel  utilization  constant.  Further,  the  outlet  fuel  temperature 
can  be  controlled  by  varying  the  air  ratio,  i.e.  the  supply  of  air 
for  cooling.  FI  and  F2  in  Fig.  2  are  static  mappings,  which  can 
be  obtained  from  the  definitions  of  fuel  utilization  and  air  ratio, 


n: 


fuel 


* 


J 


SOFC 


j 


T 


\7 


K(x,  xk)  =  exp 


(9) 


Fig.  2.  Control  scheme  of  SOFC. 
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respectively, 


A^fuel  =  FI  (J)  = 


A/2F 


(44h4+^h2+4W 


in  A/4F 

N^r  =  F2(/,  AR)  =  — ^ — J  •  AR, 

xo2 


(10) 

(ID 


where  Uf  is  the  operating  fuel  utilization.  To  meet  the  remaining 
control  objective,  what’s  needed  is  just  to  keep  the  air  ratio  in  a 
certain  range.  Therefore,  the  air  ratio  is  taken  as  a  manipulated 
variable  and  the  outlet  fuel  temperature  as  a  controlled  variable 
in  the  control  scheme.  Correspondingly  in  our  dynamic  temper¬ 
ature  model,  the  average  current  density  and  the  air  ratio  are 
two  input  variables,  and  the  outlet  fuel  temperature  is  an  output 
variable.  The  inlet  gas  temperature  and  composition  herein  are 
assumed  to  be  fixed. 

As  the  NARX  model  is  a  general  model  structure  suitable 
for  representing  a  wide  class  of  nonlinear  dynamic  systems,  it 
is  used  to  describe  the  nonlinear  temperature  dynamics  of  an 
SOFC  in  this  study.  The  dynamic  temperature  model  thereby 
can  be  written  as: 


Tk  —  fiTk—  1,  •  •  •  ?  Tk-ny,  Jk—h  •  •  •  >  Jk—nd »  AR^_i,  .  .  .  , 

AR  k-nu),  (12) 

where  Tk,  Jk,  and  AR^  represent  the  outlet  fuel  temperature,  the 
average  current  density  and  the  air  ratio  at  sampling  instant  k , 
respectively,  ny ,  rid  and  nu  are  the  maximum  lags  of  T,  J  and  AR, 
respectively,  which  are  determined  by  the  specific  process  data, 
/(•)  is  an  unknown  nonlinear  mapping.  Let  yk  =  Tk  and  Xk  = 
IX*- 1,  •  •  •  ,  Tk-ny,  Jk- 1,  •  •  •  ,  Jk-nd,  AR*_i,  .  .  .  ,  ARk-nu]T • 
Thus,  the  dynamic  temperature  modeling  problem  could  be 
stated  as:  build  an  LS-SVM  regression  model  in  the  form  of 
Eq.  (8)  based  on  a  training  set  {xk,  yk)k=\  to  approximate  the 
nonlinear  mapping /(•). 


3.3.  Construction  of  training  data  set 

Constructing  a  training  set  containing  sufficient  representa¬ 
tive  data  points  is  an  important  step  in  building  an  LS-SVM 
model  that  can  represent  SOFC’s  temperature  dynamics  over 
a  relatively  wide  operating  range.  The  basic  guideline  on  con¬ 
structing  a  training  set  is  that  the  data  in  it  should  cover  the  entire 
expected  operating  range  of  the  LS-SVM  model.  As  analyzed 
in  Section  3.2,  the  dynamic  temperature  model  has  two  input 
variables  J,  AR  and  one  output  variable  T.  Hence  the  expected 
operating  ranges  of  all  the  three  variables  should  be  covered. 

In  this  study,  the  training  data  is  generated  with  a  physical 
SOFC  model  (see  Appendix  A).  A  scheme  consisting  of  a  set 
of  operating  cases  is  designed  to  collect  training  data  by  the 
physical  model  (refer  to  Section  4.1).  Each  operating  case  is 
determined  by  three  parameters:  /,  AR  and  the  initial  tempera¬ 
ture  To.  Several  levels  of  J  and  AR,  which  distribute  uniformly 
over  their  respective  operating  ranges,  are  selected,  and  each 
level  of  J  and  AR  should  be  included  in  the  scheme.  Moreover, 
the  To  of  each  case  should  ensure  that  its  corresponding  temper¬ 
ature  dynamic  process  can  cover  the  expected  operating  range 


of  T  as  widely  as  possible.  Thus,  the  lower  or  upper  bound  of 
the  expected  operating  range  of  T is  usually  selected  as  To.  After 
resampled,  the  dynamic  data  of  these  operating  cases  can  be  used 
to  construct  the  training  set.  To  reduce  the  amount  of  training 
data,  at  first  only  the  data  of  part  of  these  cases  are  added  to  the 
training  set.  Then,  the  data  of  other  cases  are  added  gradually 
until  satisfactory  modeling  results  are  achieved.  This  approach 
guarantees  that  the  training  set  can  contain  enough  representative 
data  points  in  a  relatively  small  size. 

3.4.  Tuning  hyperparameters  automatically  hy  GA 

With  the  RBF  kernel  chosen  a  priori,  there  are  two  parameters 
that  need  to  be  tuned,  i.e.  the  regularization  parameter  y  in  Eq.  (5) 
and  the  kernel  parameter  a2  in  Eq.  (9).  The  two  parameters  are 
usually  referred  to  as  hyperparameters.  Choosing  an  optimal  set 
of  hyperparameters  is  also  an  important  step  in  LS-SVM  design. 
The  aim  of  tuning  the  hyperparameters  is  to  make  the  LS-SVM 
model  achieve  better  generalization  ability  (i.e.  the  ability  to  cor¬ 
rectly  predict  samples  that  are  not  contained  in  the  training  set), 
which  is  usually  evaluated  using  an  estimated  generalization 
error  [23].  Consequently,  the  problem  of  hyperparameter  tuning 
is  equal  to  looking  for  a  set  of  hyperparameters  in  the  search 
space  that  can  minimize  the  estimated  generalization  error.  In 
this  paper,  the  genetic  algorithm  is  employed  to  implement  the 
task  of  optima  searching.  The  simplicity  and  high  efficiency  of 
GA  can  help  to  accurately  find  the  optimal  hyperparameters  in 
a  relatively  low  time  cost. 

GA  is  a  form  of  evolutionary  algorithm,  which  is  invented  by 
John  Holland  and  has  achieved  notable  success  on  various  hard 
optimization  problems  [19,20].  The  basic  idea  of  GA  originates 
from  the  mechanisms  of  evolution  and  the  principle  of  natu¬ 
ral  selection  in  nature.  GA  starts  with  an  initial  set  of  random 
solutions  called  population.  Each  individual  in  the  population 
represents  a  feasible  solution  to  the  problem  at  hand.  The  popula¬ 
tion  “evolves”  through  successive  iterations,  called  generations. 
At  each  generation,  the  individuals  are  evaluated  using  some 
measures  of  fitness,  and  a  series  of  genetic  operations  including 
selection,  crossover  and  mutation  are  imposed  upon  individuals 
in  the  current  population  to  create  the  next  generation.  After 
several  generations,  the  algorithm  may  converge  to  the  best 
individual,  which  hopefully  represents  the  optimal  solution  to 
the  problem.  During  the  implementation  of  GA,  what’s  needed 
is  just  the  information  of  fitness.  This  is  a  prominent  differ¬ 
ence  between  GA  and  classical,  derivative-based,  optimization 
algorithms. 

For  the  problem  of  hyperparameter  tuning  by  GA,  each  set 
of  y  and  a2  is  taken  as  an  individual  in  a  population,  and  the 
estimated  generalization  error  as  the  fitness.  As  the  &-fold  cross- 
validation  is  a  very  reliable  method  to  estimate  the  generalization 
error  [23],  it  is  employed  in  this  paper.  In  k-fold  cross-validation, 
the  training  data  is  randomly  split  into  k  roughly  equal  sub¬ 
sets.  An  LS-SVM  decision  rule  is  trained  using  (k  —  1)  of  these 
subsets  and  validated  on  the  subset  left  out.  This  procedure  is 
repeated  k  times  with  each  of  the  k  subsets  used  as  the  valida¬ 
tion  subset  in  turn.  Averaging  the  validation  errors  over  the  k 
trials  gives  an  estimate  of  the  generalization  error.  The  flow- 
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Fig.  3.  Flowsheet  of  hyperparameter  tuning. 


sheet  of  tuning  the  hyperparameters  by  GA  is  illustrated  in 
Fig.  3. 

4.  Modeling  results  and  model  validation 

4.1.  Modeling  results 

A  planar  SOFC  in  a  co-flow  configuration  is  adopted  in  our 
modeling.  This  is  a  single  square  electrolyte- supported  SOFC 
operating  with  direct  internal  steam  reforming  of  methane  and 
air.  The  parameters  and  design  operating  conditions  of  the  SOFC 
are  listed  in  Table  2.  The  expected  operating  ranges  of  the  LS- 
SVM  temperature  model  are  Je  [2000,  3500]  Am-2,  AR e  [5, 
8]  and  T  e  [1247,  1347]  K,  which  include  the  design  point  of 
the  SOFC.  To  collect  training  data,  four  levels  of  J  and  AR 
are  selected,  respectively,  and  a  scheme  containing  16  operating 
cases  is  designed  (see  Table  3,  the  cases  that  are  not  used  have 
been  omitted).  By  solving  the  physical  model  using  the  fourth- 
order  Runge-Kutta  method  [27],  the  input  and  output  data  of 
the  16  cases  are  collected  and  then  resampled  using  a  sampling 
interval  of  220  s.  The  maximum  lags  ny,  nd  and  nu  in  Eq.  (12) 


Table  2 


Parameters  and  assumed  operating  conditions  of  the  SOFC 


Cell  parameters  [24-26] 

Cell  number 

i 

Cell  active  area  (m2) 

0.1  x  0.1 

Anode  thickness  (m) 

5  x  1CT5 

Cathode  thickness  (m) 

5  x  1(T5 

Electrolyte  thickness  (m) 

1.5  x  10“4 

Anode  volume  (m3) 

5.4  x  1(T6 

Cathode  volume  (m3) 

5.4  x  10-6 

Mass  (kg) 

0.2668 

Density  (kgm-3) 

6600 

Heat  capacity  (Jkg-1  K-1) 

400 

Anode  electrical  conductivity  (Q~l  m-1) 

(9.5  x  107/7)exp(- 1150/7) 

Cathode  electrical  conductivity  (£2-1  m-1) 

(4.2  x  107/7)exp(- 1200/7) 

Electrolyte  electrical  conductivity  (£2-1  m-1) 

3.34  x  104exp(-10, 300/7) 

Limiting  current  density  (Am-2) 

8000 

Activation  energy  of  anode  (J  mol-1) 

1  x  105 

Activation  energy  of  cathode  (J  mol-1) 

1.17  x  105 

Pre-exponential  factor  of  anode  (Am-2) 

5.5  x  108 

Pre-exponential  factor  of  cathode  (Am-2) 

7.0  x  10s 

Assumed  operating  conditions 

Average  current  density  (Am-2) 

3000 

Fuel  utilization 

0.85 

Air  ratio 

7 

Operating  pressure  (bar) 

1 

Inlet  gas  temperature  (K) 

1123 

Inlet  gas  composition 

Fuel:  H2,  0.2626;  CH4, 
0.171;  H20,  0.4924;  CO, 
0.0194;  C02,  0.0436 

Air:  02,  0.21;  N2,  0.79 

are  determined  by  a  trial  method,  and  herein  [ ny ,  nd,nu]T  =  [  1,1, 
1]T  is  employed.  Finally  the  training  set  {xk,  yk)l= i  C  R3  x  R 
is  obtained. 

The  two-dimensional  search  space  of  y  and  o2  is  [1, 
104]  x  [0.1,  103].  The  population  size  and  the  maximum  gener¬ 
ation  number  are  set  to  15  and  25,  respectively.  By  the  proposed 
GA-based  tuning  method  with  the  5 -fold  cross-validation  error 
as  fitness,  it  takes  7.42  h  to  find  the  optimal  set  of  (y,  a2)  at 
(8709,  33.01)  on  an  AMD  Sempron™  1.6  GHz  computer  with 
512  MB  RAM.  The  best  and  mean  fitness  values  in  each  genera¬ 
tion  are  illustrated  in  Fig.  4.  Then  an  LS-SVM  model  is  obtained 
using  the  optimal  hyperparameter  values. 


Table  3 

Operating  cases  for  constructing  the  training  set 


Case 

J  (Am”2) 

AR 

To(  K) 

Case 

/(Am'2) 

AR 

To(  K) 

1 

2000 

5 

1247 

9 

3000 

5 

1247 

2 

2000 

6 

1247 

10 

3000 

6 

1247 

3 

2000 

7 

1247 

11 

3000 

7 

1247 

4 

2000 

8 

1247 

12 

3000 

8 

1247 

5 

2500 

5 

1247 

13 

3500 

5 

1247 

6 

2500 

6 

1247 

14 

3500 

6 

1247 

7 

2500 

7 

1247 

15 

3500 

7 

1247 

8 

2500 

8 

1247 

16 

3500 

8 

1247 
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Fig.  4.  Best  and  mean  fitness  values  in  each  generation. 


4.2.  Model  validation 

First  it  is  necessary  to  note  that  there  are  two  different 
modes  of  operation  of  a  trained  LS-SVM  NARX  model:  the 
series-parallel  mode  and  the  parallel  mode  [28].  Eq.  (12)  repre¬ 
sents  the  series-parallel  mode.  The  LS-SVM  model  operating  in 
this  mode  requires  both  input  and  output  data  from  the  process 
as  inputs  and  can  predict  the  process  output  only  one  time  step 
into  the  future.  On  the  other  hand,  the  LS-SVM  model  operating 
in  the  parallel  mode  has  the  following  form: 

Tk  =  f(Tk- 1,  •  •  •  ,  Tk-ny ,  *4-1,  .  •  .  ,  Jk-nd,  AR&-1,  .  .  .  , 

ARk-nu ),  (13) 

where  the  output  of  the  LS-SVM  model  Tk  is  delayed  and  fed 
back  to  the  model  inputs  to  predict  the  process  output  further 
into  the  future.  In  this  parallel  mode,  the  LS-SVM  model  can 
be  used  independently  from  the  process  and  provide  long  range 
predictions  from  only  process  input  information.  Obviously  an 
LS-SVM  model  that  can  operate  accurately  in  the  parallel  mode 


is  a  more  powerful  tool  for  designing  control  strategies  than  that 
can  operate  accurately  only  in  the  series-parallel  mode.  How¬ 
ever,  the  parallel  mode  has  higher  requirements  for  the  prediction 
accuracy  of  the  LS-SVM  model,  because  the  prediction  error  in 
one  step  may  accumulate  to  a  bigger  one  after  multiple  steps.  To 
fully  verify  its  ability  to  represent  the  temperature  dynamics  of 
the  SOFC,  in  this  section  the  LS-SVM  model  is  operated  in  both 
the  two  modes  and  the  associated  mean  squared  errors  (MSEs) 
in  both  modes  are  calculated. 

Two  sets  of  validation  tests  (see  Table  4),  which  are,  respec¬ 
tively,  called  Validation  A  and  Validation  B  in  this  paper,  are 
carried  out  to  validate  the  LS-SVM  model.  Validation  A  is  to 
compare  the  temperature  dynamics  of  the  LS-SVM  model  with 
that  of  the  physical  model  under  operating  cases  that  are  not 
used  during  training.  The  twelve  validation  cases  in  Validation 
A  cover  the  main  operating  region  of  the  LS-SVM  model.  In 
Validation  B,  the  LS-SVM  model  is  tested  on  two  different  step 
inputs  ( AR  and  J)  and  the  corresponding  responses  are  compared 
with  those  of  the  physical  model. 

The  MSEs  of  each  validation  case  are  listed  in  Table  4.  It  is 
found  that  the  LS-SVM  model  operating  in  the  series-parallel 
mode  can  reproduce  the  temperature  dynamics  of  the  SOFC  with 
high  accuracy  for  each  validation  case.  As  expected,  the  MSE  in 
the  parallel  mode  is  bigger  than  that  in  the  series-parallel  mode 
for  each  case.  The  biggest  MSE  is  1.7255  for  Case  A12,  which 
indicates  that  the  maximum  root  mean  squared  error  (RMSE)  of 
the  temperature  prediction  is  less  than  1 .5  K.  This  is  still  very  sat¬ 
isfactory  for  multi-step  ahead  prediction.  Figs.  5-7  illustrate  the 
temperature  responses  of  the  LS-SVM  model  operating  in  the 
parallel  mode  along  with  the  comparison  with  those  of  the  physi¬ 
cal  model.  The  corresponding  figures  of  the  series-parallel  mode 
are  omitted  here  due  to  similarity  and  their  higher  accuracy.  It 
can  be  seen  from  these  figures  that  the  temperature  responses 
of  the  LS-SVM  model  show  good  consistency  with  those  of  the 
physical  model.  All  these  validation  results  demonstrate  that  the 
LS-SVM  model  is  capable  of  emulating  the  temperature  dynam¬ 
ics  of  the  SOFC  with  good  accuracy  from  only  process  input 
data. 


Table  4 

Validation  cases  and  associated  MSEs 


Case 

J  (A  m-2) 

AR 

To(  K) 

MSE  in  series-parallel  mode  ( K 2) 

MSE  in  parallel  mode  ( K 2) 

Al 

3300 

6.5 

1247 

0.0067 

0.3034 

A2 

2800 

6.5 

1247 

0.0009 

0.0327 

A3 

2300 

6.6 

1247 

0.0069 

0.4979 

A4 

3300 

6.5 

1347 

0.0097 

0.5304 

A5 

2800 

6.5 

1347 

0.0015 

0.1028 

A6 

2300 

6.5 

1347 

0.0060 

0.5439 

Validation  A 

Al 

3300 

7.5 

1247 

0.0004 

0.0085 

A8 

2800 

7.5 

1247 

0.0042 

0.1877 

A9 

2300 

7.5 

1247 

0.0067 

0.3062 

A10 

3300 

7.5 

1347 

0.0006 

0.0044 

All 

2800 

7.5 

1347 

0.0061 

0.2614 

A12 

2300 

7.5 

1347 

0.0294 

1.7255 

Case 

Input  signal  description 

MSE  in  series-parallel  mode  ( K 2) 

MSE  in  parallel  mode  ( K 2) 

Validation  B 

B1 

B2 

Step  input  of  AR  with  fixed  J= 
Step  input  of  J  with  fixed  AR  = 

3000  Am”2 

7 

0.0084 

0.0093 

0.8152 

0.8024 

Outlet  fuel  temperature  (K) 
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Fig.  5.  Temperature  responses  of  the  LS-SVM  model  operating  in  a  parallel 
Cases  A7-A9;  (d)  Cases  A10-A12. 


and  the  physical  model  for  Validation  A:  (a)  Cases  A1-A3;  (b)  Cases  A4-A6;  (c) 


Fig.  6.  Temperature  responses  of  the  LS-SVM  model  operating  in  a  parallel  Fig.  7.  Temperature  responses  of  the  LS-SVM  model  operating  in  a  parallel 

mode  and  the  physical  model  for  Case  Bl:  (a)  temperature  responses;  (b)  step  mode  and  the  physical  model  for  Case  B2:  (a)  temperature  responses;  (b)  step 

signal  of  AR.  signal  of  J. 
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5.  Conclusions 

A  dynamic  temperature  model  of  an  SOFC  is  developed  in 
this  paper  using  least  squares  support  vector  machines  (LS- 
SVMs).  Issues  concerning  the  development  of  the  LS-SVM 
temperature  model,  including  variable  selection,  training  set 
construction  and  hyperparameter  tuning,  are  discussed.  Com¬ 
prehensive  validation  tests  have  shown  that  the  developed 
LS-SVM  model  is  sufficiently  accurate  to  emulate  the  tempera¬ 
ture  response  of  the  SOFC  over  a  relatively  wide  operating  range 
using  only  process  input  information.  The  powerful  ability  of 
the  LS-SVM  temperature  model  benefits  from  the  approaches 
of  constructing  the  training  set  and  tuning  the  hyperparameters 
automatically  by  GA,  besides  the  modeling  method  itself.  In 
this  study,  the  training  data  is  generated  with  a  physical  model. 
If  sufficient  data  from  a  real  plant  is  obtained,  the  optimal  LS- 
SVM  temperature  model  of  SOFC  can  be  built  automatically  in 
the  same  way. 

The  LS-SVM  temperature  model  presented  in  this  paper  is 
useful  for  designing  temperature  control  strategies.  Control  of 
the  SOFC  using  the  LS-SVM  model  will  be  the  focus  of  our 
future  work. 
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where  kv  =  105  Pa  bar-1,  Xi  is  the  molar  fraction  of  species  /,  aq 
the  stoichiometric  coefficient  of  species  i  in  reaction  j,  and  rj  the 
rate  of  reaction  j. 

The  reaction  rates  of  (l)-(3)  are  calculated  by  the  following 
equations  [25]: 


(A.3) 


r2  =  AYkYpanX CH4  exp 


(A.4) 


r3  =  Ash&shPanCxCO*H20  —  ^sh*C02*H2),  (A.5) 


where  the  activation  energy  of  reforming  reaction 
E^r  =  82000  J  mol-1,  the  equilibrium  constant  of  WGS 
reaction  Ks^  is  obtained  using  the  method  of  Ref.  [32]. 

According  to  the  continuity  equations,  the  outlet  total  molar 
flow  rate  is  calculated  by 

Kf  =  Nfn  +  2r2,  N°*  =  NT-rj.  (A.6) 
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A.2.  Energy  balance  equation 

The  energy  conservation  law  is  used  for  the  whole  SOFC  and 
the  energy  balance  equation  [25,31]  is  written  as 


Appendix  A.  A  physical  model  of  DIR-SOFC 

This  model  is  developed  based  on  conservation  laws  and 
electrochemical  principles.  Several  assumptions  are  made  as 
follows: 


(1)  Three  reactions  (1)— (3)  are  included  in  the  model. 

(2)  The  temperature  of  the  solid  part  (including  the  PEN  and 
interconnect)  and  the  molar  fraction  of  each  gas  species 
distribute  uniformly  over  the  SOFC,  i.e.  the  SOFC  is  treated 
as  a  well-stirred  reactor  [29,30]. 

(3)  The  outlet  stream  temperature  is  equal  to  the  solid  tempera¬ 
ture:  7°ut  =  V  [29,30].  For  simplicity,  they  are  also  denoted 
by  T  in  this  paper. 

(4)  The  rate  of  energy  accumulation  in  the  gas  phase  is  negli¬ 
gible  with  respect  to  the  solid  part  [29-31]. 


A.l.  Mass  balance  equations 


,  ,  dr  A  A 
mscs„ — +yy 
p  dr 


k=li=l 


d rtihi 
d t 


Q  gen? 


(A.7) 


where  ft  =  5  for  CH4,  H2,  CO,  C02,  H20,  ft  =  2  for  02,  N2. 

The  first  term  on  the  left  side  of  Eq.  (A.7)  is  the  rate  of  energy 
accumulation  in  the  solid  part  of  the  SOFC  and  the  second  term, 
which  can  be  neglected  because  of  the  fourth  assumption,  is 
that  in  the  gases.  The  first  term  on  the  right  side  is  the  enthalpy 
change  due  to  the  input  and  output  gas  streams,  the  second  term 
is  the  heat  generated  by  the  three  reactions  and  is  calculated  by 

2  gen  =  AH\ri  +  +  A/73T3  —  Pfc.  (A.8) 

The  gas  specific  heat  cpi  in  Eq.  (A.7)  is  a  function  of  temperature: 

cpi{T)  =  b0  +  hT  +  b2T2  +  b3T\  (A.9) 


The  molar  balance  is  applied  to  each  gas  species  within  the  ^  3  Voltage  model 
SOFC  as  follows: 


Anode  : 


Pan  Vin  dX[ 

RT  dT 


j= 1 


(A.l) 


The  output  voltage  is  calculated  by  subtracting  three  kinds 
of  overpotential  losses,  i.e.  ohmic,  activation  and  concentration 
overpotentials,  from  the  reversible  voltage  as  Eq.  (A.  10). 

U  =  C/n  —  7? ohm  —  7? act  —  ^7conc 


i  e  {CH4,  H2,  CO,  C02,  H20}, 


(A.  10) 
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The  reversible  voltage  is  calculated  by  the  Nernst  equa¬ 
tion: 

RT 

UN  =  U0(T)  +  —  In 

with  Uo(T)  =  1.2723-2.7645  x  10~4T,  where  Uo(T)  is  the 
reversible  voltage  at  standard  pressure  [33]. 

The  ohmic  overpotential  is  calculated  by  the  Ohm’s  law  as: 

7? ohm  =  JRohm,  (A.  12) 


0.5  0.5 
XK2X02  Pea. 

xh2o 


(A.ll) 


where  R0hm  is  the  internal  resistance  of  the  cell  and  is  calculated 
from  the  conductivity  of  the  individual  layers  [6]  by 


D  Tin  Tlec  Ta 

^ohm  — - 1 - 1 - • 

^an  ^elec  ^ca 

The  activation  overpotential  is  given  by: 


(A. 13) 


Oact 


2  RT 
neF 


sinh 


-l 


2j0,a 


+ 


2  RT 
neF 


sinh 


-l 


2jo,c 


(A.  14) 


where  jo, an  and  jo,Ca  denote  the  anodic  and  cathodic  exchange 
current  density,  respectively,  and  are  modeled  in  the  following 
form  [34]: 

(^)  exp  (-|f)  ,  (A- 15) 

» <»  =  ft* ( )  “p(  If  )  <A-i6) 

The  concentration  overpotential  is  determined  from 
RT  f  J\ 

0  conc  =  -— lnfl--H,  (A.17) 

where  /l  is  the  limiting  current  density. 
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